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ABSTRACT 

The HIFI instrument on board the Herschel Space Observatory has been used to observe interstellar nitrogen hydrides along the sight-line towards 
G 10.6-0.4 in order to improve our understanding of the interstellar chemistry of nitrogen. We report observations of absorption in NH N=l <— 0, 
7 = 2 «— 1 and ortho-NH2 L,i «— Oo,o- We also observed ortho-NH^ lo «— Oo, and 2o «— lo, para-NH^ 2\ <— li, and searched unsuccessfully for 
NH + . All detections show emission and absorption associated directly with the hot-core source itself as well as absorption by foreground material 
over a wide range of velocities. All spectra show similar, non-saturated, absorption features, which we attribute to diffuse molecular gas. Total 
column densities over the velocity range 1 1 - 54 km s _I are estimated. The similar profiles suggest fairly uniform abundances relative to hydrogen, 
approximately 6 x 10~ 9 , 3 x 10~ 9 , and 3 x 10~ 9 for NH, NH2, and NH3, respectively. These abundances are discussed with reference to models of 
gas-phase and surface chemistry. 

Key words. ISM: abundances - ISM: molecules - Submillimetre: ISM - Molecular processes - Line: formation - Astrochemistry 
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1. Introduction 

Molecular hydrides are important in the chemistry of the inter- 
stellar medium since they often appear in the first steps in chains 
of reactions that lead to other more complex species. The pro- 
duction pathways of nitrogen-bearing molecules are still rather 
uncertain since key species, such as NH + , NH, and NH2 have not 
been widely observed. Even the first identified polyatomic in- 
terstellar molecule, ammonia (NH3), has been widely observed 
mainly in its para symmetry form, which leaves its formation 
mechanism poorly constrained in the diffuse molecular gas. Both 
gas-phase chemistry and grain surface reactions have been pro- 
posed as formation mechanisms, but clearly more observations 
are needed for a better understanding. 

Both NH and NH2 are well known in comets (e.g. Swings et 
al. 1941; Meier et al. 1998; Feldman et al. 1993), and have been 
observed in stellar photospheres (e.g. Schmitt 1969; Farmer & 
Norton 1989) via their electronic, vibration-rotation, and high 
rotational transitions. Interstellar NH was first detected in the in- 
terstellar medium by Meyer & Roth (1991) by optical absorption 
spectroscopy. Subsequent observations by Crawford & Williams 
(1997) and Weselak et al. (2009) have yielded six lines of sight 
where column densities of both NH and H2 are directly mea- 
sured. The average value of the column density ratio in these 
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* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



diffuse and translucent sightlines is Af(NH)/W(H 2 ) = 3 x 10~ 9 . 
Interstellar NH2 was first observed by van Dishoeck et al. ( 1 993) 
in absorption towards Sgr B2 in three fine-structure components 
of the para-NH2 li,o - lo,i transition with partially resolved 
hyperfine structure at frequencies 461 to 469 GHz. Further ab- 
sorption lines of both ortho and para forms of NH2 and NH 
were observed through use of the long-wavelength spectrome- 
ter aboard the Infrared Space Observatory (ISO: Cernicharo et 
al. 2000; Goicoechea et al. 2004; Polehampton et al. 2007). The 
ISO observations were unable to resolve the hyperfine structure 
of either molecule. Ammonia, NH3, was the first polyatomic 
molecule to be identified in interstellar space (Cheung et al. 
1968) by means of its microwave inversion transitions. Although 
NH3 has been widely observed in dark clouds and star-forming 
regions, there have been very few measurements of the inversion 
lines in the diffuse interstellar gas (Nash 1990; Liszt et al. 2006). 

With the unique capabilities of Herschel Space Observatory 
(Pilbratt et al. 2010) transitions between 157 and 625 yum (0.48- 
1 .9 THz) are available with the Heterodyne Instrument for the 
Far-Infrared (HIFI; de Graauw et al. 2010) with very high sensi- 
tivity. This allows searches for spectrally resolved, ground-state 
rotational transitions of NH + , NH, NH2, and NH3 with the same 
instrument. Abundances of these species are key diagnostics for 
the nitrogen chemistry. 

The most sensitive and model-independent method for 
measuring column densities of interstellar molecules is high- 
resolution absorption line spectroscopy, which is being ex- 
ploited by the Herschel PRISMAS key programme (PRobing 
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InterStellar Molecules with Absorption line Studies). This pa- 
per presents the first observations and analysis of nitrogen hy- 
drides along the sight-line towards the massive star-forming re- 
gion G10.6-0.4 (W31C) as a part of the PRISMAS program. 
Observations of interstellar absorption towards seven additional 
submillimetre-wave continuum sources will be presented and 
analysed in later publications. 

The ultra-compact HII region G10.6-0.4 in the star-forming 
W31 complex is an extremely luminous submillimetre and in- 
frared continuum source. The source is located within the so- 
called 30 kms -1 arm at a kinematic distance of 4.8 kpc (Fish 
et al. 2003). The gas associated directly with G10. 6-0.4 is de- 
tected in the velocity range Vlsr ~ -15 to +5kms~ 1 , while the 
foreground gas along the line of sight is detected at Vlsr ~ 5 to 
55kms _1 . The focus of this paper is on the diffuse interstellar 
gas traced by the absorption lines. 



2. Observations 

The observations, which took place in March 2010, are sum- 
marised in TableQ] The 2q-1o and 2i-li NH3 transitions were 
observed in the same band. We used the dual beam switch 
mode and the wide band spectrometer (WBS) with a band- 
width of 4 GHz and an effective spectral resolution of 1.1 MHz. 
The corresponding velocity resolution is about 0.3 kms -1 at the 
higher frequencies and 0.6 km s 1 at 572 GHz. In addition, si- 
multaneous observations were performed using the high- resolu- 
tion spectrometer (HRS) with an effective spectral resolution of 
0.25 MHz (Av ~ 0.1 km s" 1 ) and a bandwidth of 230 MHz. Three 
observations were carried out with different frequency settings 
of the local oscillator (LO) corresponding to a change of approx- 
imately 15kms~' to determine the sideband origin of the lines. 
Two orthogonal polarisations were used during all observations. 



Table 1. Observed transitions. 
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"The NH and NH 2 frequencies refer to the strongest hyperfine struc- 
ture component. *HIFI consists of 7 different mixer bands and two 
double sideband spectrometers. All transitions were observed in the 
upper sideband except NH + . r System temperature. rf The on-source 
time for each integration. 



The half-power beam width of the telescope is 37, 22, and 
17" at 572, 953, and 1215 GHz, respectively. Pointing was cen- 
tered at RA. = 18 h 10 m 28 s :70, Dec. = -19? 55! 50" (72000). The 
reference beams were located 3' on either side of the source and 
the calibration is described by Roelfsema (2010). 

We have used the standard Herschel pipeline to Level 2 
which provides fully calibrated spectra. The data were first anal- 




Velocity [km s ] 



Fig. 1. Double sideband spectra of NH, 0-NH2, o and P-NH3, 

and NH + over the LSR velocity range -50 to 85 km s~' . 



ysed using Herschel interactive processing environment (HIPEq 
Ott 2010) version 2.4, and in parallel we also used the software 



1 HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia. 
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package xfl The data quality is excellent with very low inten- 
sity ripples, two polarisations that are in agreement to better than 
10%, and very good agreement between the three LO-tunings 
without contamination from the image sideband. We thus aver- 
age the three LO-tunings and both polarisations in both spectra, 
except for NH3 at 572 GHz and NH + which display ripples in 
the V- and H-polarisation, respectively, and were therefore not 
included. The resulting rms noise is 11, 120, 50, 51, and 74 mK 
for NH 3 at 572 and 1215 GHz, for NH, NH 2 , and NH + , respec- 
tively. 



3. Results 

Figure[T|shows the double sideband WBS spectra of all observed 
transitions. The frequency scale is converted as customary to 
Doppler velocities relative to the local standard of rest (Vlsr)- 
All spectra, except NH + which is not detected, show emission at 
negative velocities and a deep absorption at -1 km s associated 
with G10.6-0.4 as measured by observations of OH masers. At 
higher velocities we see a number of absorption features which 
are also seen in atomic hydrogen (Fish et al. 2003), HCO + , HNC, 
HCN and CN (Godard et al. 2010). The only transition which 
does not show absorption components from the foreground ma- 
terial is ortho-NHi 2q- 1q. 

The double-sideband calibrated antenna temperatures shown 
in Fig. [T]have to be divided by two in order to get the correct 
continuum levels, ^(cont), which are 0.6, 2.6. 2.6, 3.0 and 3.4K 
for NH 3 at 572 GHz, NH 2 , NH, NH + , and NH 3 at 1 215 GHz, re- 
spectively. The sideband gain ratio has been shown to be close to 
unity in previous PRISMAS observations (Neufeld et al. 2010b; 
Gerin et al. 2010), although departures at the 10% level have 
been seen at some frequencies (Neufeld et al. 2010a). We here 
adopt a sideband ratio of unity. 

Figure [2] shows a comparison of NH, NH2, and para- 
NH3 where the intensities are normalised to the continuum in 
a single sideband as rA/7A(cont)-l assuming a sideband gain 
ratio of unity. The NH and NH2 spectra are strikingly simi- 
lar, despite their complicated hyperfine structures. The para- 
NH3 spectra also show the same absorption pattern, although 
not as strongly. Figure [3] shows a similar comparison of ortho- 
and para-NHi,. Here, too, we see the same absorption pattern 
in both species, and strong emission and self-absorption in the 
orf/10-line from the source itself. The strongest velocity compo- 
nents in the foreground material lie at Vlsr = 16, 22.5, 27.5, 30 
and 40.5 km s" 1 . 

The line opacities are calculated as 

T=-ln(7'ii ne /rA(cont) - 1). To be conservative and not in- 
clude absorption from the source itself, we use the velocity 
interval ll-54kms~ 1 . The total integrated opacities are 
/ rdV = 20, 16, 9, and 3.4km s" 1 for NH, NH 2 , NH 3 at 572 and 
1215 GHz, respectively. 

Even though unsaturated absorption lines provide straight- 
forward determinations of opacity and column densities, the 
numerous hyperfine-structure components and velocity features 
complicate the analysis of both NH and NH 2 . The spectro- 
scopic properties of NH and NH 2 have been studied extensively 
and references, rest frequencies and molecular properties for 
NH and NH 2 are listed in the Cologne Database for Molecular 
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Fig. 2. Comparison of NH, o-NH 2 and p-NH 3 at 1 215.2 GHz 

where the intensities have been normalised to single sideband 
continuum. 
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Fig. 3. Comparison of o- and p-NH 3 at 572.5 and 1 215.2 GHz 
where the intensities have been normalised to single sideband 
continuum. 

Spectroscopy (CDMS, Miiller et al. 2005), and NH 3 in Jet 
Propulsion LaboratorjQ (JPL, Pickett et al. 1998). 

In order to estimate the total column density from the inte- 
grated opacity, it is necessary to correct for the population of 
molecules in unobserved excited levels and to quantify any pos- 
sible effect of stimulated emission. We have done this through 
non-equilibrium excitation calculations of integrated opacity 
with an enhanced version of the RADEX codqj (van der Tak 
et al. 2007). Where collision rates are unknown, we have made 
guesses scaled in proportion to radiative line strengths. We as- 
sume diffuse molecular cloud conditions with 7\ = 30 K and 
«(H 2 ) =500crrT 3 . The background radiation is a model of the 
average Galactic background radiation in the solar neighbour- 
hood plus the cosmic microwave background radiation. The ex- 
citation of the hydride molecules in diffuse gas is dominated by 



2 Developed by Per Bergman at Onsala Space Observatory, Sweden; 
http : //www . Chalmers . se/rss/oso-en/ observations/ 



3 http://www.cdms.de 

4 http://spec.jpl.nasa.gov 
5 



The published version will soon include these enhancements, see 
http : //www . sron . rug . nl/~vdtak/radex/index . shtml 
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the background continuum radiation and not by collisions. The 
resulting excitation temperatures of the observed submillimetre 
transitions are typically 4-5 K, which are small enough com- 
pared to hv/k that no correction for emission is required. All 
the optical depths of the hyperfine components that contribute 
to the line are summed. We have treated ortho- and para-forms 
of molecules together in the expectation that their relative abun- 
dances are fixed by the chemical formation process. Plausible 
rates of formation and destruction are included explicitly in the 
excitation calculation. Thus we derive total column densities as 
summarised in Table [2] The analysis is complicated for NH3 by 
the existence of metastable excited states, especially in the para 
species. The entry for NH3 in Table [2] is the sum of the mea- 
sured column density in Jk - Oo, the measured column density 
in the lower inversion sublevel of 1 1, and calculated values in the 
upper inversion sublevel of li and in both sublevels of 22. All 
other states are expected to have negligible populations. The un- 
certainties of measurement and in the excitation analysis do not 
fully constrain the ortho/ para ratio in NH3, but it appears that 
the effective formation temperature must be greater than 20 K, 
otherwise the strength of the 1215 GHz line cannot be explained. 

Table 2. Total column densities, N, and abundances with respect 
to hydrogen, X. 



Species x 


N" 
(cm" 2 ) 


X = NJN B 


NH 


1.5xl0 14 


5.6xl0- 5 


NH 2 


8.0xl0' 3 


3.0x10"" 


NH 3 


8.7x10" 


3.2x10-* 


NH + 


<2xl0 13 


<8xl0- 10 



"Calculated in the velocity interval 1 1 - 54 km s 



We note that the three neutral nitrogen hydrides have com- 
parable abundances: NH is approximately twice as abundant as 
NH2 and NH3. To estimate the abundance with respect to hy- 
drogen we use Nh =2.7xl0 22 ctrT 2 as total column density of 
hydrogen in all forms (Neufeld et al. 2010b). This value is con- 
sistent with the estimated atomic hydrogen column density in 
this line-of-sight, N(H) = 1.2xl0 22 cirT 2 (Godard et al. 2010) 
inferred from absorption profiles observed with the VLA (Fish 
et al. 2003), and implies that the absorbing material is mostly 
molecular even if there most likely are large variations between 
the individual velocity components. Our estimate of the abun- 
dance of NH relative to total hydrogen is somewhat higher 
than the NH/H2 ratio of 3 x 10~ 9 for six lines of sight where 
both are observed directly. The resulting abundances are found 
in Table [2] The non-detection of NH + is used to put an up- 
per limit on its column density and abundance (3<x). A previ- 
ous limit on NH + towards f Oph from ultraviolet observations, 
N(NH + )/Nh < 1 x 10~ 6 , is orders of magnitude less sensitive 
than our limit (Snow 1979; de Almeida & Singh 1982). 

4. Discussion 

The production of NH and NH2 by purely gas-phase processes 
is problematic since their formation is inhibited by a lack of a 
sufficient source of N + . The atomic ion can be formed by cos- 
mic ray ionization of N or by reaction of He + with N2 or CN, 



which were formed by neutral-neutral reactions. The standard 
gas-phase ion-molecule chemistry of nitrogen hydrides is then 
initiated by the reaction N + +H2 — » NH + + H, because H^ does 
not react rapidly with N. Typical steady state dark cloud mod- 
els («=lxl0 3 -5xl0 4 cnr 3 , T= 10-40 K, A v > 10), predict 
an NH2 abundance of (1-I0)xl0~ 8 , an NH abundance 10 times 
lower, and an NH2/NH3 ratio of 0.3 - 1 .5 for a wide range of as- 
sumptions (e.g. Millar et al. 1991; Langer & Graedel 1989), but 
these are not directly applicable to diffuse molecular gas. 

Previous observations of NH in diffuse clouds were taken to 
imply that NH is produced on grain surfaces (Mann & Williams 
1984; Wagenblast et al. 1993; O'Neill et al. 2002). Such mod- 
els, however, often predict up to 1 000 times more NH3 than 
NH 2 (Hasegawa & Herbst 1993). An additional source of NH 
and NH2 is fragmentation of ammonia by photodissociation 
in diffuse clouds. Models of ultraviolet illuminated PDRs pre- 
dict NH 2 /NH 3 <1 and NH 3 /NH<1 in regions with extinction 
A v < 5 mag (e.g. Sternberg & Dalgarno 1995). The abundance 
patterns that we see in diffuse molecular gas are clearly different 
from those in Sgr B2, where NH 3 /NH 2 /NH ~ 100/10/1 and the 
fractional abundance of NH is a few times xlO ~ 9 (Cernicharo 
et al. 2000; Goicoechea et al. 2004). 

We have used two different approaches to chemical mod- 
els in order to compare the nitrogen-hydride abundances with 
theory. First, we constructed a pseudo-time-dependent model 
with constant physical conditions including both gas-phase and 
surface chemistry (Hasegawa et al. 1992) with Ay -2 and 3, 
n H = 500 and 5 000 cm 3 , T gas = 30 and 50 K, and r dust = 10 K. 
The ultraviolet field is 1 Go, and the cosmic ray ionization 
rate is g= 1.3xlO" 17 S -1 . The resulting abundances are found 



in Fig. IA.1I (online material). Gas-phase chemistry alone does 
not match the observed abundances, failing by factors of 10 
to 100 for NH and NH2. Inclusion of processes on dust sur- 
faces increases the abundances of the three neutral hydrides to 
10- 8 -10~ 9 , but fails to match the high NH/NH3 ratio. If the 
dominant source of NH3 is association on surfaces and if the 
ortho/ para ratio is fixed at a surface temperature below 20 K, 
then the 1215 GHz para line should be much weaker relative to 
the ortho 572 GHz line than is observed. The observed line ratio 
is perhaps somewhat better explained by a high formation tem- 
perature, which points to exoergic gas-phase formation of NH3. 
A second approach is to construct models of weak photon- 
dominated regions (PDR) through use of the Meudon PDR 
code (Le Petit el al. 2006; Goicoechea & Le Bourlot 2007), 
which is a steady-state, depth-dependent model with pure 
gas-phase chemistry. We use parameters appropriate for dif- 
fuse clouds, UV = 0.5-5 Go, n B = 100-1 OOOcirT 3 , f = 10~ 17 - 
10~ 15 s , Ay < 5 and illuminate one side only. The results of 



the best fitting model are found in Fig. IA.2I (online material). 
Although the resulting relative abundances of nitrogen-hydrides 
approach the observed ratios inside the cloud, the abundances 
relative to hydrogen are still too low. 

5. Conclusions 

Our first detections of spectrally resolved rotational transitions 
of nitrogen hydrides in the interstellar medium show remarkable 
similarities of line profiles and abundances towards one back- 
ground source. The formation mechanisms are, however, dif- 
ficult to explain. Neither standard gas-phase nor grain-surface 
chemistry can fully explain our observations, which may suggest 
that important reactions have been overlooked in the chemistry 
of interstellar nitrogen. Further analysis of the line profiles will 
yield abundance ratios in separate velocity components and will 
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allow the nitrogen species to be compared directly with tracers 
of molecular hydrogen like CH and HE Our sample will also 
be enlarged with seven additional sight-lines in the PRISMAS 
project. The chemical models will be refined and the possible 
role of turbulent dissipation regions will be assessed. 
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Fig. A.l. A pseudo-time-dependent model with constant physical conditions taking into account both gas-phase and grain surface 
chemistry. For all lines: Tjust = 10 K and the solid/dashed lines show r gas = 30 and 50 K, respectively. The black thin lines show the 
observed abundances, and the dashed black line the NH + upper limit. Black line: gas phase & grain surface chemistry using Ay - 3, 
«h = 5000 cirT 3 . Red line: gas phase & grain surface chemistry using Ay -2, tin- 500cm -3 . Grey line: only gas phase chemistry 
using Ay — 3, «h = 5000 cirT 3 . Blue line: only gas phase chemistry using Ay = 2, tin- 500cm -3 . 
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Fig. A.2. The cloud depth-dependent Meudon PDR model for pure gas-phase chemistry using «h = 1 000 cm , £= 10" - 10" s 
and a UV-field of 10 times the interstellar radiation field in Draine units (~5 Go). The left panel shows the relative abundances of 
NH/NH2 and NH2/NH3, and the right panel show the absolute abundances. 



